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Abstract 

This article addresses the steady three-dimensional flow of an Oldroyd-B nanofluid over a bidirectional stretching surface 
with heat generation/absorption effects. Suitable similarity transformations are employed to reduce the governing partial 
differential equations into coupled nonlinear ordinary differential equations. These nonlinear ordinary differential equations 
are then solved analytically by using the homotpy analysis method (HAM). Graphically results are presented and discussed 
for various parameters, namely, Deborah numbers /}, and fi 2 , heat generation/absorption parameter X, Prandtl parameter 
Pr, Brownian motion parameters^, thermophoresis parameter A', and Lewis number Le. We have seen that the increasing 
values of the Brownian motion parameter N, and thermophoresis parameter N, leads to an increase in the temperature 
field and thermal boundary layer thickness while the opposite behavior is observed for concentration field and 
concentration boundary layer thickness. To see the validity of the present work, the numerical results are compared with the 
analytical solutions obtained by Homotopy analysis method and noted an excellent agreement for the limiting cases. 
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Introduction 

During the past few years, study of the boundary layer flow of 
nanofluids over a linear stretching surface has become more and 
more attractive because of its numerous applications in industrial 
manufacturing. With regard to the sundry application of 
nanofluids, the researchers have been given considerable attention 
to improve heat transfer using nanofluids. Regular fluids, such as 
ethylene, water, glycol mixture and some types of oil have low heat 
transfer rates. Therefore it is necessary to improve some physical 
properties such as thermal conductivity and heat transfer rate of 
conventional fluids by the utilization of nanoparticles in base fluid. 
The term nanofluid was first time introduced by Choi [1]. In 
another paper, Choi et al. [2] observed that thermal conductivity 
of pure fluid can be increased by a factor of 2 with an addition of 
one percent by volume fraction of the nanoparticle. 

Sakiadis [3] was the first who investigated the boundary layer 
flow on a continuous stretching surface. In his paper, he provided 
numerical solutions of the boundary layer flow over a continuous 
stretching surface. Later on Crane [4] analyzed the exact solution 
of boundary layer flow of Newtonian fluid due to stretching of an 
elastic sheet moving linearly in its own plane. Wang [5] 
investigated the free convection on a vertical stretching surface. 
Heat transfer analysis over an exponentially stretching continuous 
surface was analyzed by Elbashbeshy [6]. Rana and Kango [7] 
discussed the effect of rotation on thermal instability of compress- 
ible Walters' (model) elastico-viscous fluid in porous medium. Heat 
transfer over a stretching surface with variable heat flux in 



micropolar fluids was presented by Ishak et al. [8] . Chamkha and 
Aly [9] examined MHD free convective boundary layer flow of a 
nanofluid along a permeable isothermal vertical plate in the 
presence of heat source or sink. Thermosolutal convection in 
Walters' (Model B') elastico-viscous rotating fluid permeated with 
suspended particles and variable gravity field in porous medium in 
hydromagnetics was investigated by Rana [10]. Matin et al. [11] 
presented the MHD mixed convective flow of a nanofluid over a 
stretching sheet. Chand and and Rana [12] examined the 
oscillating convection of nanofluid in porous medium. Aziz and 
Khan [13] studied natural convective flow of a nanofluid over a 
convectively heated vertical plate. Kuznetsov and Nield [14] 
analyzed the natural convective flow of a nanofluid past a vertical 
plate. Khan and Pop [15,16] investigated the laminar flow of 
nanofluids past a stretching sheet. Hamad et al. [17] formulated 
the problem of free convective flow of nanofluid past a semi- 
infinite vertical plate with influence of magnetic field. Hady et al. 
[18] investigated the effects of thermal radiation on the viscous 
flow of a nanofluid and heat transfer over a non-linear sheet. 
Makinde and Aziz [19] performed the numerical study of 
boundary layer over a linear stretching sheet. Cheng [20] analyzed 
the behavior of boundary layer flow over a horizontal cylinder of 
elliptic cross section in a porous. Narayana and Sibanda [21] 
elaborated the effects of laminar flow of a nanofluid over an 
unsteady stretching sheet. Kameswaran et al. [22] investigated 
flow due to a stretching or shrinking sheet with viscous dissipation 
and chemical reaction effects. The effects of an unsteady 
boundary-layer flow and heat transfer of a nanofluid over a 
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porous stretching/shrinking sheet have been investigated by 
Bachok et al. [23]. Hamad and Ferdows [24] presented the 
similarity solutions for viscous flow and heat transfer of a nanofluid 
over a non-linear stretching sheet. The studies on heat generation/ 
absorption effects for boundary layer flow of nanofluids are very 
limited. Recently, Alsaedi et al. [25] investigated the effects of heat 
generation/ absorption on stagnation point flow of nanofluid over 
a surface with convective boundary conditions. Thermal instability 
of Rivlin-Ericksen Elastico-Viscous nanofluid saturated by a 
porous medium were presented by Chand and Rana [26]. On 
the onset of thermal convection in rotating nanofluid layer 
saturating a Darcy-Brinkman porous medium were studied by 
Chand and Rana [27]. Nandy and Mahapatra [28] examined the 
effects of slip and heat generation/ absorption on MHD stagnation 
point flow of nanofluid past a stretching/ shrinking surface with 
convective boundary conditions. On the onset of thermosolutal 
instability in a layer of an Elastico-Viscous nanofluid in porous 
medium was investigated by Rana et al. [29]. 

However, to the best of author's knowledge, no attempts have 
thus far been communicated with regards to free convective 
boundary layer flow of three-dimensional Oldroyd-B nanofluid 
over a stretching surface. The aim of the present article is to study 
the free convective boundary-layer flow of three-dimensional 
Oldroyd-B nanofluid fluid flow over a stretching sheet. The 
Oldroyd-B fluid model was employed to describe rheological 
behavior of viscoelastic nanofluid. The Oldroyd-B fluid model is 
important because of its applications in the production of plastic 
sheet and extrusion of polymers through through a slit die in 
polymer industry. The considered stretched flow problem involves 
problem involves the significant heat transfer between the sheet 
and the surrounding fluid. The extrudate in this mechanism starts 
to solidify as soon as it exits from the die and then sheet is collected 
by a wind-up roll upon solidification. Physical properties of the 
cooling medium, e.g., its thermal conductivity has pivotal role in 
such process. The success of whole operation closely depends upon 
the viscoelastic character of fluid above the sheet. By applying 
boundary layer approximations a system of nonlinear partial 
differential equations is obtained. Then, invoking suitable similar- 
ity transformations, we reduced the system into nonlinear ordinary 
differential equations. This system of coupled nonlinear ordinary 
differential equations is then solved analytically by using the 
homotpoy analysis method (HAM). The variations of different 
flow controlling parameters on the velocity, temperature and 
concentration profiles are addressed. 

Mathematical Formulation 

Consider a steady three-dimensional (x,y, z) free convective 
boundary layer flow of an incompressible Oldroyd-B nanofluid 
over a stretching sheet kept at a constant temperature T w and 
concentration C w . The ambient temperature and concentration 
far away from the sheet are taken as T x and C x , respectively. 
The flow is due to a bidirectional stretched surface at z = 0. The 
governing equations for the steady three-dimensional flow of an 
Oldroyd-B nanofluid, approximated by boundary-layer theory, 
are [30] 
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Here (u, v, w) are velocity components, T the temperature, C 
the concentration, Ai and the relaxation and retardation times 
respectively, p the fluid density, a the thermal diffusivity, Qo the 
heat generation/absorption parameter, x the ratio of effective heat 
capacity of the nanoparticle material to the heat capacity of the 
fluid, Db the Brownian diffusion coefficient and Dt the 
thermophoresis diffusion coefficient. 

Equations (1) to (5) are subjected to the following boundary 
conditions 



u = ax, v = by,w = 0,T= T w , C = C w at z = 0, 



(6) 
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Table 1. Convergence of the homotopy solutions when p A = /? 2 =0.2, /? = 0.4, Pr =1.2,JVj=JV t = 0.1 andL e = l are fixed. 





Order of approximation 


-f "(0) 


-g"(0) 


-8(0) 


-<p'(0) 


1 


1 .00420 


0.337840 


0.622000 


0.352000 


5 


1.02196 


0.328912 


0.549080 


0.493576 


10 


1.02155 


0.328848 


0.549423 


0.489147 


15 


1.02154 


0.328870 


0.549446 


0.488882 


20 


1.02154 


0.328869 


0.549438 


0.488934 


26 


1.02154 


0.328869 


0.549438 


0.488939 


30 


1.02154 


0.328869 


0.549438 


0.488939 


35 


1.02154 


0.328869 


0.549438 


0.488939 
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f"' + (f+g)f"-f' 2 + ^ [2(f+g)f'f"-(f+g) 2 f") + 
lh[2(f"+g")f"-(f+g)f]=0, 

g"' + (f+g)g"-g 2 + lh W+g)g'g"-(J' + g) 2 g"'] + 
P 2 [2(f"+g")g"-(f+g)g ir ]=0, 

9" + Pr (f+g)ff+ PrjVif 9'+ Pr N t 9' 2 + Prl0 = O, (11) 



«-»0, v-»0, T->T X , C->C X as z->oo. (7) 

The similarity variables are introduced as 
u{x,y) = axf'(r\), v(x,y)=ayg'(}i), w= -y/av\f(ri) + g(r])} 

and Eqs.(l)-(7) can be cast as 




Figure 1. Variation of /i, on %) when p 2 =0.2, p = 0A, Pr =0.4, 2 = 0.2, N h =N, = 0.1 and Le = 1.0 are fixed. 

doi:1 0.1 371 /journal.pone.01 051 07.g001 
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Figure 3. Variation of /? on %) when /?, = £,=0.2,/? = 0.4, Pr =1.2, X = 0.2, N b =N, = 0.1 and Le=1.0 are fixed. 
doi:1 0.1 371 /journal.pone.01 051 07.g003 



PLOS ONE | www.plosone.org 



4 



August 2014 | Volume 9 | Issue 8 | e105107 



Oldroyd-B Nanofluid 




0 2 4 6 8 

Figure 4. Variation of Pr on 0(jj) when ft = /J 2 =0.2, /J = 0.2, X = 0.2,N l =N, = 0.1 and Le=1.0 are fixed. 
doi:10.1371/journal.pone.0105107.g004 
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Figure 5. Variation of X on %) when /}, =/J 2 = 0.2, /i = 0.4, Pr = 1.2, jV 4 =JV, = 0.1 and Le=1.0 are fixed. 
doi:10.1371/journal.pone.0105107.g005 
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Figure 7. Variation of N,, on 0(n) when ft = ft = 0.2, /J = 0.4, Pr =1.2, X = 0.2, N, = 0.1 and Le=1.0 are fixed. 
doi:10.1371/journal.pone.0105107.g007 
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Figure 8. Variation of N, on %) when ft = ft = 0.2,/? = 0.4, Pr = 1.2, X = 0.2,N h = 0.1 and Le=1.0 are fixed. 

doi:10.1371/journal.pone.0105107.g008 




0 2 4 6 8 



Figure 9. Variation of ft on %) when ft = 0.2,/? = 0.4, Pr = 1.2, 1 = 0.2, JV 4 = AT, = 0.1 and£e=1.0 are fixed. 
doi:10.1371/journal.pone.0105107.g009 
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Figure 10. Variation of ft, on %) when ft =0.2,/? =0.4, Pr = \.2,X = 0.2,N b =N, = Q.\ and Le = 1.0 are fixed. 

doi:10.1371/journal.pone.0105107.g010 
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Figure 11. Variation of X on ,9(ij) when ft = ft = 0.2,0 = 0.4, Pr = 1.2,% =JV« = 0.1 and Le=1.0 are fixed. 

doi:10.1371/journal.pone.0105107.g011 
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Figure 13. Variation of N b on Sfa) when ft =0 2 = 0.2, 0 = 0.4, Pr = 1.2,^ = 0.2,^ = 0.1 and Le = 1.0 are fixed. 

doi:10.1371/journal.pone.0105107.g013 
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Figure 14. Variation of N, on %) when 0, =#j = 0.2,/? = 0.4, Pr = 1.2, 1=0.2, Ni, = 0.1 and Le=1.0 are fixed. 
doi:10.1371/journal.pone.0105107.g014 



Pre(f+gW+^6" = 0, 



f=0,g=0,f=l,t 



: 1,5=1 at »7 = 0, 



(12) 



(13) 



/'-►0,g'->0, e->o,^-»o as r/- 



(14) 



where prime denotes differentiation with respect to r\. Moreover, 
Pi and /?2 are the Deborah numbers, /? the ratio of stretching rates 
parameter, Pr the generalized Prandtl number, X the heat source 
().>0) and the heat sink (a<0) parameter, Nb the local Brownian 
motion parameter, N r the local thermophoresis parameter and he 
the Lewis number which are defined as 



Table 2. A comparison for the velocity gradients for different values of [1 when ji x 



= 0 are fixed. 



p 


HPM result [31] 


HPM result [31] 


Exact result [31] 


Exact result [31] 


Present result 


Present result 


-f "(0) 


-g"{0) 


-f"(0) 


-g"(0) 


-f "(0) 


-g"(0) 


0.0 


1.0 


0.0 


1.0 


0.0 


1.0 


0.0 


0.1 


1 .02025 


0.06684 


1 .020259 


0.066847 


1.02026 


0.06685 


0.2 


1 .03949 


0.14873 


1 .039495 


0.148736 


1.03949 


0.14874 


0.3 


1.05795 


0.24335 


1.05794 


0.243359 


1.05795 


0.24336 


0.4 


1.07578 


0.34920 


1.075788 


0.349208 


1.07578 


0.34921 


0.5 


1.09309 


0.46520 


1.093095 


0.465204 


1.09309 


0.46521 


0.6 


1.10994 


0.59052 


1.109946 


0.590528 


1.10994 


0.59053 


0.7 


1.12639 


0.72453 


1.126397 


0.724531 


1.12639 


0.72453 


0.8 


1.14248 


0.86668 


1.142488 


0.866682 


1.14249 


0.86668 


0.9 


1.15825 


1.01653 


1.158253 


1.016538 


1.15826 


1.016538 


1.0 


1.17372 


1.17372 


1.173720 


1.173720 


1.17372 


1.17372 
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Table 3. Comparision of results for the local Nusselt number —6'(0) in the absence of non-Newtonian parameters and 
nanoparicles when [i = 0 with the work of Khan and Pop [15] and Nadeem and Hussain [32] 





Pr 


Present result 


Khan and Pop [15] 


Nadeem and Hussain [32] 


0.07 


0.066 


0.066 


0.066 


0.20 


0.169 


0.169 


0.169 


0.70 


0.454 


0.454 


0.454 


2.0 


0.911 


0.911 


0.911 
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The physical quantities of interest are the local Nusselt number 
Nu x and the local Sherwood number Sh x , which are defined as 



Nu x -- 



Sh x -- 



8T 

5z 



dC 



(16) 



(17) 



In terms of dimensionless form one has 

Re iNu x = -&{G), ReiSh x =- ${<!>), (18) 
where Re = ux/v is the local Reynolds number. 

Convergence of the Homotopy Solutions 

The problems containing non-linear coupled ordinary differen- 
tial equations (19)— (12) subjected to boundary conditions (13) — (14) 
have been computed analytically by the homotopy analysis 
method (HAM). In the HAM role of the auxiliary parameters ly, 
fig, fig, and h v is of key importance because they control the 
convergence of the series solution. The most suitable value of these 
auxiliary parameters is calculated by considering minimum square 
error which is given by 



N + l 



E 

,/=o 



N f J2Fj(iAri) 



(19) 



Table 1 ensure the convergence of the series solution which 
shows that the convergent solution for the velocity is obtained at 
20th-order of approximation whereas such a convergence for 
temperature and concentration is achieved at 26th-order of 
approximation. 

Numerical Results and Discussion 

The aim of this section is to analyze the influence of the various 
physical parameters on the velocity, temperature and nanoparticle 
fields respectively. Figs. 1-14 are plotted to see the variation of the 
Deborah numbers ft and ft, Prandtl number Pr, heat source 
(1>0) or sink (a<0) parameter, Brownian motion parameter Nt 
and thermophoresis parameter N t on the fluid temperature and 
concentration fields. 

Fig. 1 shows the influence of the Deborah number ft on the 
temperature field. By increasing Deborah number ft both the 
fluid temperature and thermal boundary layer thickness increases. 
This is due to fact that the Deborah number ft involves relaxation 
time X\ . An increase in the relaxation time leads to increase in the 
temperature and boundary layer thickness. Fig. 2 illustrates the 
effects of the Deborah number ft on the temperature field. From 
this figure, it is noted that the behavior of the Deborah number ft 
is opposite to that of ft. This is due to fact that the retardation 
time provides resistance which causes a reduction in the 
temperature and thermal boundary layer thickness. Fig. 3 presents 
the effects of the stretching parameter ft on the fluid temperature 
d(r]). We observed that the temperature and thermal boundary 
layer thickness reduce with the increasing ft Fig. 4 illustrates the 
influence of the Prandtl number Pr on the temperature field. We 



Table 4. Comparision of results for the local Nusselt number -9'(0) and local Sherwood number -<j>'(0) in the presence of 
nanoparticle of when /? = 0, ft = ft = 0.3, /l = 0, Pr =6 and Le= 1 are fixed with the work of Nadeem et al. [33]. 





N, 


N b 


Present result 


Present result 


Nadeem et al. [33] 


Nadeem et al. [33] 


-9(0) 


f'(0) 


-9(0) 


<|>'<0) 


0.3 


0.3 


0.33984 


1.83994 


0.33988 


1.83935 


0.5 


0.3 


0.24088 


1.95813 


0.24099 


1.95862 


0.3 


0.5 


0.14814 


1 .87029 


0.14820 


1.87035 


0.5 


0.5 


0.10478 


1 .94565 


0.10486 


1.94572 
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observed that the temperature and thermal boundary layer 
thickness are reduced for large Prandtl number. Since thermal 
diffusivity is an agent which plays a key role for lower or higher 
temperature. Hence resulting larger value of the Prandtl number 
corresponds to diminishing of the thermal diffusivity resulting in a 
temperature decrease. Figs. 5 and 6 are plotted to analyze the 
effects of the heat source parameter ( when X>0) and heat sink 
parameter ( when X<0), respectively. It is observed that 
temperature of the fluid increases with the increase in the heat 
source parameter and an opposite behavior is observed for heat 
sink parameter. The influence of the Brownian motion parameter 
Nb and thermophoresis parameter N, on the temperature is 
depicted through Figs. 7 and 8, respectively. It is observed that the 
temperature and thermal boundary layer thickness increases as the 
Brownian motion parameter Nb increases. Physically, this is due to 
the fact that with an increase of the Brownian motion parameter 
Nb the random motion of particle increases which results in an 
enhancement in the temperature profile. The temperature and 
thermal boundary layer thickness are detected to increase with an 
increase in thermophoresis parameter N, (Fig. 8). In fact with the 
increase of the thermophoresis parameter N, the difference 
between the wall temperature and reference temperature increases 
which causes increase in temperature profile. 

In Figs. 9 and 10, we plotted the concentration profile for 
various values of the Deborah numbers f5\ and P 2 , respectively. As 
the Deborah numbers Pi increases, the concentration profile as 
well as concentration boundary layer thickness increase. However, 
the effects of p 2 on the concentration profile are quite opposite to 
that of /?j . Fig. 1 1 shows the influence of the heat generation 
parameter X on the concentration profile. A decrease in the 
concentration profile and concentration boundary layer thickness 
near the plate is noted while the reverse effect is reported far away 
from the plate with the increasing value of the heat generation 
parameter X. Fig. 12 illustrates the influence of the Lewis number 
he on the concentration profile (j>(ri). It is noted that the 
concentration profile increases by increasing the Lewis number 
he as Lewis number is inversely proportional to the diffusion 
coefficient. Thus an increase in Lewis number yields a decrease in 
diffusion which finally results in a decrease of mass fraction 
function 9(rj). The variations with r\ of the concentration profile 
for different values of the Brownian motion parameter Nb and 
thermophoresis parameter Nt are presented in Figs. 13 and 14, 
respectively. In Fig. 13, it is observed that concentration profile 
increases with the increasing of the Brownian motion parameter 
Nb. This is due to the dependency of the concentration on the 
temperature field and we expect that a lower Brownian motion 
parameter allow a deeper penetration of the concentration. On the 
other hand, a qualitatively opposite trend in the concentration 
profile is observed as the thermophoresis parameter N, increases. 
Further, it is noticed that the thermophoresis parameter A^ affects 
the concentration profile more than Brownian motion parameter 
Ni, does. 

Numerical values for the velocity gradients —f"(0) and —g"(0) 
axe compared with the existing literature in the absence of both 
nanoparticles and non-Newtonian effects and shown in table 2, 
where they are found to be in excellent agreement, cementing the 
validity of the present results. Table 3 gives comparison of local 
Nusselt number — (?'(0) with the results obtained by Khan and Pop 
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the involved parameters. It is reported that the local Nusselt 
number — 6'(0) increases when P and Pr increase whereas it 
decreases as X,Nb, N t and he increase. It is evident from table 5 
that the local Sherwood number —<j)'(0) increases with the 
increase of the parameters p, Pr , X, N[, and he, however, it 
decreases with the increase of N t . 

Concluding Remarks 

This study has analyzed the effects of the heat generation/ 
absorption on three-dimensional flow of an Oldroyed-B nanofluid 
over a bidirectional stretching sheet. From the present investiga- 
tion, the main observations were as follows: 

• Qualitatively, effects of the Deborah numbers P\ and Pi on the 
temperature and concentration profiles were similar. 

• The temperature profile as well as thermal boundary layer 
thickness were increased by increasing both the Brownian 
motion parameter Ni, and thermophoresis parameter N t . 

• The temperature of the fluid and thermal boundary layer 
thickness is enhanced when there is a increase in the heat 
generation parameter X. 

• The concentration profile was decreased with the increase of 
the Brownian motion parameter A^ and a quite opposite 
behavior was noted with increasing thermophoresis parameter 
N t . 

• It was noted that the thermophoresis parameter N, affected the 
concentration profile more than the Brownian motion 
parameter Nb did. 

• An increase in the heat generation parameter (whenl>0) 
corresponds to reduction in the values of the local Nusselt 
number — 9'(0) while the opposite behavior is observed for the 
local Sherwood number — (j)'(0). 

• The magnitude of the local the local Nusselt number — 8'(0) 
decreases with the increase of the Brownian motion parameter 
N b . 

• The magnitude of the local the local Sherwood number — (j)'(0) 
increase with the increase of the Brownian motion parameter 
N h . 
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